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ABSTRACT 

In the present work we use a deep-exposure spectrum taken by the SUMER spectrometer in a polar coronal 
hole in 1996 to measure the ion temperatures of a large number of ions at many different heights above the 
limb between 0.03 and 0.17 solar radii. We find that the measured ion temperatures are almost always larger 
than the electron temperatures and exhibit a non-monotonic dependence on the charge-to-mass ratio. We use 
these measurements to provide empirical constraints to a theoretical model of ion heating and acceleration 
based on gradually replenished ion-cyclotron waves. We compare the wave power required to heat the ions 
to the observed levels to a prediction based on a model of anisotropic magnetohydrodynamic turbulence. We 
find that the empirical heating model and the turbulent cascade model agree with one another, and explain the 
measured ion temperatures, for charge-to-mass ratios smaller than about 0.25. However, ions with charge-to- 
mass ratios exceeding 0.25 disagree with the model; the wave power they require to be heated to the measured 
ion temperatures shows an increase with charge-to-mass ratio (i.e., with increasing frequency) that cannot be 
explained by a traditional cascade model. We discuss possible additional processes that might be responsible 
for the inferred surplus of wave power 

Subject headings: line: profiles — Sun: corona — Sun: UV radiation — techniques: spectroscopic — turbu- 
lence — waves 



1. INTRODUCTION 

In order to understand the physical processes that heat the 
solar corona and accelerate the solar wind, theories must be 
able to predict measurements of the plasma parameters in the 
regions that are being heated and accelerated. In the low- 
density, open-field regions that correspond to polar coronal 
holes on the Sun, the plasma is expected to become colli- 
sionless very near the solar surface, and thus the individual 
particle species (e.g., protons, electrons, and heavy ions) may 
exhibit diverging properties. These differences are key probes 
of the microscopic kinetic physics that drives the heating and 
acceleration. The Solar Ultraviolet Measurements of Emit- 
ted Radiation (SUMER) spectrometer aboard the Solar and 
Heliospheric Observatory iSOHO) has measured these prop- 
erties for more than a decade (see, e.g., WiUielm et al. 1995, 
1997). 

Inspired by both the SUMER observations presented below 
and other remote-sensing and in situ measurements, we fo- 
cus on the damping of ion cyclotron resonant Alfven waves 
as a natural means of preferentially heating and accelerating 
coronal ions. Efficient kinetic interactions between ion cy- 
clotron waves and particle velocity distributions in the solar 
wind have been studied for several decades (e.g., Abraham- 
Shrauner & Feldman 1977; Hollweg & Turner 1978; Marsch 
et al. 1982; Hollweg 1986; Tu & Marsch 1997; Cranmer et 
al. 1999a; Hollweg & Isenberg 2002; Cranmer 2002; Kohl 
et al. 2006). Despite much intense theoretical work, the ori- 
gin of the high-frequency (10^-10"* Hz) cyclotron waves in 
the solar corona is not yet known. Most ideas for the genera- 
tion of these waves fall into two classes: (1) excitation at the 
base of the corona, and (2) gradual growth or replenishment 
of cyclotron waves throughout the corona and solar wind. The 
first idea has been proposed to occur from impulsive processes 



taking place in magnetically complex regions (e.g., reconnec- 
tion in microflares; see Axford & McKenzie 1992; Ruzmaikin 
& Berger 1998). The second scenario has been proposed to 
occur mainly from magnetohydrodynamic (MHD) turbulent 
cascade (Isenberg & Hollweg 1983; Cranmer & van Balle- 
gooijen 2003; Chandran 2005), but other processes such as 
kinetic microinstabilities have also been suggested. Indirect 
support for the second idea of "gradual replenishment" has 
come from empirical studies of radio scintillations (Hollweg 
2000) and theoretical models of the charge and mass depen- 
dence of the wave damping (Cranmer 2000, 2001). In our 
paper, we will base our model on the gradual replenishment 
scenario. 

One of the main predictions of ion-cyclotron models is pref- 
erential heating and acceleration of ions in the solar wind, 
with the effect of waves on ions being dependent on the Z,/A, 
ratio, where Z, is the ion's charge and A, its atomic mass (with 
both usually specified in proton units). One of the observables 
that can best constrain theoretical models of ion-cyclotron 
waves is the behavior of the ion temperature (i.e., the mag- 
nitude of ion random motions on microscopic scales). These 
models are also dependent on the overall amplitude of MHD 
waves in the corona, which can also be constrained empir- 
ically by the so-called non-thermal component of emission 
line profiles. 

Measuring ion temperatures and non-thermal velocities is 
however a very difficult task, as these are two separate phys- 
ical quantities that have to be determined from one single 
observable, the line width. Thus, some assumptions have 
to be made. However, these limit the accuracy of the mea- 
surements. In the past, most assumptions were taken on the 
value of Ti, since the quantity of interest was the non-thermal 
velocity Vnth- The most common assumptions were Ti = Tg, 
or Ti = Tm, where Tm is the temperature of formation of the 
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line, usually taken as the temperature of maximum abun- 
dance of the ion emitting the line. A more elaborate set of 
assumptions were adopted by Seely et al. (1997), who as- 
sumed that ions formed in the same temperature range have 
the same 7^; this allowed them to determine both 7/ and Vnth 
from SUMER/50//0 observations. Tu et al. (1998), on the 
contrary, assumed that Vnth is the same for all ions. This al- 
lowed them to set an upper limit to v„,i„ and conservatively 
assumed Vnth = as a lower limit. From these limits, the range 
of possible 7} values could be determined. In doing this, Tu 
et al. (1998) limited the number of assumptions to the mini- 
mum and chose not to determine Vnth, so that they could set 
rather precise limits to T, of the ions in their dataset. Dolla & 
Solomon (2008) went a step further by making assumptions 
not on the values of 7} or Vnth, but on the physical process that 
affects them. They assumed that: ( 1 ) non-thermal velocities 
are due to undamped Alfven waves, and (2) that Mg X was not 
heated by them so that its 7] was constant and could be used 
to measure Vnth- 

The aim of the present work is twofold. First, we will ex- 
tend the analysis done by Tu et al. (1998) by applying their 
method to spectra emitted by coronal hole plasmas at the north 
pole, that include many more lines from many more ions than 
considered by them. We will repeat the 7) measurements for 
as many heights above the solar limb as made possible by 
the signal-to-noise ratio and instrument scattered light. This 
has been achieved by using deep-exposure SUMER spectra 
encompassing the entire wavelength range (500-1500 A) in- 
cluded by detector B on SUMER. 

Second, we will use the measured ion temperature values 
as constraints for semi-empirical theoretical models of ion 
heating and acceleration based on the damping of ion cy- 
clotron resonant Alfven waves. We will develop a model that 
is mainly intended to be illustrative of the proper ranges of 
plasma parameters that are consistent with the observations in 
the low corona, but not yet self-consistent enough to be con- 
sidered a definitive or unique explanation of the SUMER data 
set.' 

The first part of the paper deals with the ion temperature 
measurements. Section |2] introduces the method of analysis 
and Section|3]describes the observations studied in the present 
work. Our measurements are reported in Section H] The sec- 
ond part of the paper describes the comparison with the theo- 
retical model. In Section|5]the separation of line widths into 
thermal and non-thermal parts is discussed in more detail. The 
physical processes assumed in the semi-empirical models are 
given in Section|6] and a comparison with a purely theoretical 
model of MHD turbulence is made in Section [T] The model 
results are shown and described in Section |8] Section |9] sum- 
marizes the present work. 

2. METHOD OF ANALYSIS 

In this work, we have applied the analysis method described 
by Tu et al. (1998), that we briefly summarize here. The 
measured full width half maximum (FWHM) d\ of a spec- 
tral line can be related to the dynamical status of the emitting 
ion through the equation 
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' Specifically, these models do not try to validate or disprove the assump- 
tion that the ion cyclotron waves are created by an anisotropic turbulent cas- 
cade; they only try to work out the implications of such an assumption. 



where d\ is the FWHM of the line after the instrumental 
broadening has been removed, 7} and M,- = WpA,- {nip is the 
proton mass) are the ion's temperature and mass, is Boltz- 
mann's constant, c is the speed of light, Ao is the rest wave- 
length of the line, and Vnth is the ion's non-thermal velocity. 
From equation ([T]), the ion temperature given by each line can 
not be larger than the value obtained when Vnth = 0: 
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where A, is the atomic mass. In the same way, the non-thermal 
velocity can not be larger than the value obtained when 7^ = 0: 
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In order to determine a lower limit to 7], it is necessary to 
make an assumption on Vnth- The method devised by Tu et al. 
(1998) allows us to determine the minimum Ti values compat- 
ible with a set of measured FWHM values by assuming that 
Vnth is the same for all ions. If we define Vmin as the small- 
est among the Vmax./ values obtained from a set of lines, the 
ion temperature of each ion must be larger than the value ob- 
tained when v'nth = v'min- If we indicate with /' the line for which 
Vnth = Vmin, the range of ion temperatures compatible with the 
FWHM measured with line / is given by 
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Note that the lower limit of Ti for line j is zero. The inter- 
val defined by equation (HI is not to be considered the "error 
bar" due to experimental uncertainties in the measurement of 
the line's FWHM; it indicates the range of equally possible 7} 
values allowed by the observed FWHM value. Uncertainties 
in the measured dX are propagated to the limit values of 7) 
in equation (|4|, and widen the range of possible ion tempera- 
tures. This method of analysis has the advantage of reducing 
the assumptions on 7) and Vnth to the very minimum. However, 
it only provides us with an upper limit to Vnth- 

3. OBSERVATIONS AND DATA REDUCTION 

3.1. Observations 

Observations were taken on 1996 November 3, and con- 
sist of a full spectral scan encompassing the entire SUMER 
wavelength range carried out over a polar coronal hole. The 
center of the slit field of view was (0",1 149"), so that coronal 
hole plasma was observed at heliocentric distances between 
1.03 solar radii {Rq) and 1.33 Rq. The instrument configu- 
ration allows the recording of the solar spectrum in ~ 43 A- 
wide subsections of the 500-1600 A spectral range; in the 
present observation the full SUMER range was observed in 
61 subsections, each shifted from the previous one by ~ 13 A. 
This ensured that all lines were observed on the more sensitive 
KBr-coated portion of the detector at least once. Each spec- 
tral subsection was observed with SUMER detector B with 
1200 s exposure time, so that the final effective wavelength 
range of the dataset we studied is 500-1500 A. An image 
of the SUMER field of view superimposed to EIT images is 
shown in Figure [T] 
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Fig. 1 . — Field of view of the SUMER slit, superimposed to FIT images in 
all four filters: Top left: He II (5 X 10"* K); Top right: Fe IX-X (6.3 X 10^ K); 
Bottom left: Fe XII (1.4 X 10^ K); Bottom right: FeXV(2.2x 10* K). The FIT 
images were taken on the same day as SUIVIFR (3 November 1996) between 
1 UT and 2 UT. 



TABLF 1 

Lines used for the present study. Smax (in units of solar radii 

FROM DISK CENTER) INDICATES THE MAXIMUM DISTANCE BEYOND 
WHICH THE LINE COULD NOT BE USED DUE TO SCATTERED LIGHT 
CONTAMINATION. 



Ion 


Wvl. (A) 


^max 


Seq. 


log Tmax 


Zi/Ai 


NalX 


681.72 


1.171 


Li 


5.92 


0.348 


Si IX 


694.90 


1.151 


C 


6.05 


0.285 


ArVIII 


700.25 


1.100 


Na 


5.61 


0.175 


MglX 


706.06 


1.171 


Be 


5.99 


0.329 


Fe VIII 


721.26 


1.049 


K 


5.88 


0.125 


MglX 


749.55 


1.171 


Be 


5.99 


0.329 


Mg VIII 


762.65 


1.059 


B 


5.91 


0.288 


Mg VIII 


769.38 


1.079 


B 


5.91 


0.288 


Ne VIII 


770.42 


1.171 


Li 


5.80 


0.347 


Mg VIII 


772.29 


1.090 


B 


5.91 


0.288 


Ne VIII 


780.34 


1.171 


Li 


5.80 


0.347 


Mg VIII 


782.37 


1.151 


B 


5.91 


0.288 


CalX 


821.22 


1.049 


Mg 


5.80 


0.200 


Mg VII 


854.64 


1.079 


C 


5.80 


0.247 


Mg VII 


868.08 


1.110 


C 


5.80 


0.247 


Ne VII 


887.24 


1.028 


Be 


5.71 


0.297 


Ne VII 


895.16 


1.100 


Be 


5.71 


0.297 


S VI 


933.41 


1.079 


Na 


5.38 


0.156 


Ne VII 


973.33 


1.028 


Be 


5.71 


0.297 


FeX 


1028.04 


1.161 


CI 


6.08 


0.161 


Ovi 


1031.93 


1.171 


Li 


5.58 


0.313 


O VI 


1037.63 


1.171 


Li 


5.58 


0.313 


Si VII 


1049.26 


1.079 





5.79 


0.214 


Al VII 


1053.88 


1.059 


N 


5.79 


0.222 


Al VII 


1056.79 


1.059 


N 


5.79 


0.222 


Al VIII 


1057.86 


1.069 


C 


5.94 


0.259 


Cax 


557.76 


1.110 


Na 


5.87 


0.225 


Ne VI 


558.59 


1.028 


B 


5.63 


0.248 


Ne VI 


562.80 


1.028 


B 


5.63 


0.248 


Cax 


574.01 


1.110 


Na 


5.87 


0.225 


MgX 


609.80 


1.171 


Li 


6.04 


0.370 


Nv 


1238.82 


1.079 


Li 


5.24 


0.286 


FeXII 


1242.00 


1.141 


P 


6.22 


0.197 


Nv 


1242.81 


1.069 


Li 


5.24 


0.286 


MgX 


624.97 


1.161 


Li 


6.04 


0.370 


Si VIII 


1440.51 


1.028 


N 


5.92 


0.249 


Si VIII 


1445.76 


1.161 


N 


5.92 


0.249 


FeXI 


1467.07 


1.100 


S 


6.14 


0.179 



1.00 1.05 MO MS '.20 t.25 '.30 '.35 

Distance from the limb (Rg^J 

Fig. 2. — Intensity of the Mg IX 706.06 A line as a function of distance 
from the limb (in counts). The scattered Hght intensity is superimposed as a 
dashed line, and it has been normalized to match the Mg IX intensity averaged 
over the three positions farthest from the solar limb. 

3.2. Data reduction and scattered light 

Raw data were reversed and corrected from flat field and 
geometrical distortions using the standard SUMER software 
available in the SolarSoft package for IDL. Corrections for 
deadtime and local gain depression were not necessary be- 
cause count rates were always sufficiently low. A large set 
of lines was selected in order to sample as wide a range of 
Z,/A, ratios as possible with bright, isolated lines that mini- 
mized uncertainties tied to line profile fitting. These lines are 
listed in Table [T] In order to increase the signal-to-noise, the 
SUMER slit was divided into 10 pixel-wide bins along the slit 
direction: within each bin, the emission was averaged. The 
profile of each line was fitted with a Gaussian profile, obtain- 
ing line centroid, width and total intensity. 

In the present work, we are primarily interested in line 
widths, but line intensities were used to determine the 
amount of scattered light contamination. Instrument-scattered 
light consists of disk radiation scattered by mirror micro- 
roughness, that forms on the detector a ghost, unshifted spec- 
trum of the solar disk that contaminates true local emission. 
A significant scattered light blending may alter line profiles 
in ways that are not easy to assess and correct. Before at- 
tempting an analysis of line widths, it is important to check 
that the scattered light contribution is negligible and, if not 
negligible, to remove it. Scattered light can be distinguished 
from true coronal emission by considering the variations of 
the measured intensity with distance from the limb, because 
true emission decreases much faster than instrument-scattered 
light. In order to evaluate the amount of scattered light con- 
taminating each line, we first determined the rate of decrease 
of scattered light intensity with distance from the limb. This 
was done first by measuring the normalized intensity versus 
height relationship (henceforth scattered light profile) of a few 
chromospheric and transition region lines which, outside the 
solar disk, are almost entirely made of scattered light. We then 
assumed that the emission of the selected coronal lines mea- 
sured at the farthest distance from the limb is entirely com- 
posed of scattered light, and used the scattered light profile to 
evaluate the absolute intensity of instrument-scattered emis- 
sion for each coronal line at all heights. An example of com- 
parison between scattered light and true coronal emission is 
given in Figure|2]for the MglX line at 706.06 A. It is impor- 
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tant to note that this method actually overestimates the scat- 
tered light emission. We considered the scattered light to be 
negligible when its intensity, as estimated using this method, 
is 25% or less of the measured intensity of a coronal line. 

This method allowed us to determine the maximum height 
Rmax below which the scattered light can be safely considered 
negligible. For the strongest lines, we could stretch our mea- 
surements up to 1.17 Rq, while only the closest bin provided 
intensity for the weakest lines. Table[T]also reports the value 
of Rm-dx for each line. 

The presence of scattered light in line widths can also be 
checked independently by considering the quantity d\/Xo- In 
fact, from equation [T| it is easy to see that c/A/Aq only de- 
pends on physical properties of the ion and so it is the same 
for all lines of the same ion, unless some problem affects the 
profile of a line. We calculated this quantity for all the ions 
in Table [T] with more than one line in order to check whether 
there was agreement at all heights up to Rnax- In the few cases 
where dX/Xo was not constant within the same ion, we have 
investigated the cause of the disagreement, usually due to poor 
signal-to-noise or to blending lines and rejected the positions 
where there were problems. It is important to note that ions 
with only one line are considered less reliable because their 
width dX could not be checked using this method. 

4. ION TEMPERATURES 

The ion temperatures measured from line widths are shown 
in Figure |3] which displays them as a function of the Z,/A, 
ratio for a few heights above the limb. The analysis method 
provides rather wide ranges of temperature for almost all ions, 
so that any height dependence is nearly washed away and the 
results are qualitatively the same up to the maximum height 
1.17 Rq. The ions emitting weak lines only provide a few 
measurements close to the limb, and in all cases they are com- 
patible with a constant 7] value even if their 7] ranges seem to 
decrease with height. Brighter lines provide 7] measurements 
up to 1.17 Rq, that sometimes seem to increase and some- 
times seem to decrease with height; in all cases they too are 
compatible with constant T,. As height increases, fewer and 
fewer lines are useful to measure 7}, so that it is more difficult 
to sample the Z,/A,- range as accurately as with observations 
closer to the limb. 

Landi (2008) measured the electron temperature of the 
same emitting plasma under consideration in the present 
work, finding that it is nearly isothermal at each height, with 
temperature slightly increasing with height from log(7^./K) ~ 
5.90 to log(re/K) ~ 6.05. Figure|3]shows the measured upper 
and lower limits of T^, on top of the 7} ranges: in all cases the 
latter are larger than the former and if agreement is found, it is 
only because the 7} ranges are rather broad. This behavior is 
apparent even at the lowest heights and it confirms the earlier 
results by Tu et al. (1998): coronal hole plasmas are not in 
equilibrium. 

Figure [3] shows that there is a relationship between T/ and 
the Z,/A, ratio. At larger heights the number of available ions 
is limited so the relationship between 7} and Z, /A, is less ap- 
parent, or it can not even be determined, although results are 
consistent with those at lower heights. Figure |3] shows an U- 
shaped 7} versus Z,/A,- relationship, where the ions provid- 
ing large 7] at Z,/A,- < 0.2 ai-e Fex to Fexil and Arvill. 
For Z,/A, > 0.2, 7) values are more constant, showing only 
a mild increase with Z,/A,. It is interesting to note that all 
the ions with 7} closer to, or in agreement with, all have 
Zj/Aj between 0.2 and 0.3. The relationship between 7} and 
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Fig. 3. — Ion temperature versus the charge-over-mass ratio Z,/A,, for 
a selected subset of measured heights. The horizontal lines represent the 
minimum and maximum value of the electron temperature as measured by 
Landi (2008). 
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Fig. 4. — Value of Vmin as a function of height above the limb. 



Zi/Aj can be compared with the measurements carried out by 
Landi (2007) on three different quiet Sun regions using the 
same diagnostic method. Differences are striking, as quiet 
Sun Ti values show no dependence on the Z,/A, ratio, and 
they are in much better agreement with the electron temper- 
ature than in the present work. However, we found that the 
lines we selected have a sufficiently large signal-to-noise to 
provide uncertainties in the limit values of Tj that are small 
when compared to the range of 7] given by equation (|4]i. 

The Tu et al. (1998) diagnostic technique does not allow us 
to measure Vnth, since equation ^ only allows us to set an 
upper limit on it. The value of Vmin increases only slightly 
with height from ~ 30-35 km s"' to ~ 45-50 km s~', and it 
is shown in Figure |4] 

5. THERMAL AND NON-THERMAL COMPONENTS 

As discussed above, the basic assumption in the interpreta- 
tion of the off-limb line profiles is that there are only two con- 
tributors to the (mainly Gaussian) widths: "thermal" micro- 
scopic motions and "non-thermal" bulk motions along the line 
of sight (LOS) due to unresolved waves or turbulence. An- 
other key assumption is that the non-thermal part of the line 
width has no intrinsic dependence on ion charge or mass. This 
is valid when the spectrum is dominated by low-frequency 
MHD waves, but other types of dispersive waves have been 
suggested to play a role (e.g., Moran 2002; Ofman et al. 
2005). However, we retain the standard assumption of no 
charge or mass dependence of the non-thermal speed, because 
in most models of MHD turbulence the power that resides at 
kinetic dispersive scales tends to be much smaller than that at 
the larger non-dispersive scales. 

Figure |5] compares a range of observational and theoretical 
values for the non-thermal velocity, including the upper limits 
derived above. This quantity is interpreted as the transverse 
root mean square Alfven wave amplitude ((5v^) The wave 
amplitudes from theoretical models were divided by V2 in 
order to sample only the motions along one of the two per- 
pendicular directions (i.e., only along the LOS). Cranmer & 
van Ballegooijen (2005) discussed the observational data of 
Banerjee et al. (1998) and Esser et al. (1999), as well as other 
empirical constraints on Alfven wave properties. The discrep- 
ancy between the two theoretical model curves is due to the 
fact that the Cranmer & van Ballegooijen (2005) model uti- 
lized a fixed empirical set of background plasma parameters 
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Fig. 5. — Assembled non-thermal velocities, interpreted as perpendicular 
Alfven wave amplitudes and plotted versus heliocentric distance in units of 
the solar radius, from empirical non-WKB wave inodels (dot-dashed line; 
Craniner & van Ballegooijen 2005) and turbulence-driven coronal heating 
inodels (solid line; Cranmer et al. 2007). Also, empirical ineasureinents froin 
Banerjee et al. (1998) (blue crosses), Esser et al. (1999) (yellow region), and 
from this paper (red diamonds). 



(e.g., wind speed, Alfven speed, density) and Cranmer et al. 
(2007) computed these properties self-consistently with the 
waves. 

Because in coronal holes (close to the limb) the magnetic 
field appears to be mainly pointed in the radial direction, the 
LOS component of the thermal temperature (for ion /) is as- 
sumed to be identical to the perpendicular temperature r^,. 
This assumption was validated by test calculations that took 
account of the full LOS integration and used anisotropic bi- 
Maxwellian distributions. These calculations showed that the 
dominant contribution to a collisional line profile (i.e., sensi- 
tive to density squared) comes from a very narrow range of 
heights, all close enough to the plane of the sky that the sam- 
pled component of the distribution always remained close to 
Tn- 

The model calculations below mainly utilize the ion tem- 
perature data at a single reference radius. This radius was cho- 
sen as 1 .069 Rq, which seemed to be the best balance between 
having a large number of available ions (which drops off as 
height increases) and the existence of significant "preferential 
heating" away from thermal equilibrium (which is lessened 
as height decreases). The measured maximum non-thermal 
speed at this height is 34.6 km s"', which falls between the 
modeled values of about 30 km s~' (Cranmer & van Balle- 
gooijen 2005) and 40 km s"' (Cranmer et al. 2007). However, 
the actual non-thermal speed must be less than the upper limit 
derived above. Somewhat arbitrarily, we choose 25 km s"' as 
a fiducial non-thermal speed for this height, which is used to 
subtract from the total line width to obtain the individual ion 
temperatures T±i} 

^ To some extent, when comparing the ion temperatures to one another at 
the same height, the assuined value for the non-thermal wave amplitude is not 
so important. As long as the non-thermal speed is independent of ion charge 
and mass, the relative differences between the various ion temperatures are 
unchanged. 
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6. PHYSICAL PROCESSES IN THE MODEL 

The basic idea of the semi-empirical models to be discussed 
below is that we assume a turbulent replenishment of ion cy- 
clotron wave power in the corona with a known radial de- 
pendence. The absolute normalization for this wave power, 
though, is determined for a given ion of charge Z, and mass 
A, iteratively by finding the value that produces the observed 
level of perpendicular heating at the reference radius of 1 .069 
Rq. By determining this wave power for the ions, each inde- 
pendently of the others, we construct a model of the overall 
frequency dependence of the ion cyclotron power spectrum. 

The solution for the radial dependence of T±i uses an inter- 
nal energy conservation equation similar to that used by Cran- 
mer et al. (1999a) with some minor changes to be discussed 
below. The properties of the background proton-electron 
plasma are held fixed at values determined by comparison 
with empirical measurements (mainly those of Cranmer & 
van Ballegooijen 2005). The heavy ions are thus treated as 
"test particles" that do not influence the primary plasma or 
one another. A dominant effect that needs to be included is the 
temperature equilibration due to Coulomb collisions between 
protons and the heavy ion species of interest. The high den- 
sities in the low corona give rise to rapid collisions that can 
mask the preferential ion heating from ion cyclotron waves 
(by driving T±j back down to the proton temperature Tp). 

The time-steady energy conservation equation that we solve 
is given by 



dr 



TndA _ Qn + CnpiT±p-Tn) 



A dr 



niUiks 



(5) 



where the two terms on the right-hand side represent resonant 
heating iQ±i) and Coulomb collisions with protons (C±ip). 
The second term on the left-hand side accounts for adiabatic 
expansion, and depends on the radial variation in the cross- 
sectional area A of the polar flux tube. One notable depar- 
ture from the models of Cranmer et al. (1999a) is that we do 
not consider the corresponding equation for the parallel ion 
temperature T^i,-. Thus, for simplicity, we assume that there 
is no significant dependence of Tj_i on T^^j. Even for strong 
anisotropy, the isotropization (coupling) term J±ip from Cran- 
mer et al. (1999a) is a factor of ^ 0.2/A, weaker than the stan- 
dard proton-ion collision term, and thus is negligible. 

The proton-ion collision coefficient in equation ^ is given 
by C±ip = Ik^riiVip, where n,- is the ion number density. The 
Coulomb collision rate itself is 



3 



niiinpa^^a^ 



(6) 



(e.g., Spitzer 1962), and 

a|| = [{2kBT\\i/mi) + (2kiiTy/mp)y^^ 



aj_ = [(2kBT±i/mi) + (2kBT±p/mp)] 



1/2 



(7) 
(8) 

(see also Barakat & Schunk 1981, 1982; Isenberg 1984). As 
in Cranmer et al. (1999a), we take the Coulomb logarithm 
to be a constant of InA = 21. We assume the protons are 
isotropic, and we also set 7]|, = Tp in the definition of fl|| 
(which was verified to hold at the SUMER heights in trial runs 
using the original equations from Cranmer et al. 1999a). We 
do not consider collisional interactions between ion species / 
and any other heavy ion species, which implies that ion / does 
not ionize or recombine to any other charge states over the 
range of heights (see, e.g., Lie-Svendsen & Esser 2005). 



The resonant cyclotron heating rate is given in the so-called 
"optically thin" (narrow resonance zone) limit of Cranmer 
(2000) as 



niirii 



47r^r{ni,r) 



Q^Vl ( 1 



(9) 



where Bq is the radial background magnetic field strength, 
is the ion gyrofrequency, and Va is the Alfven speed. The 
resonant wave power spectrum V{oj) differs from the power 
spectrum used by Cranmer et al. (1999a) by a factor of An 
because the latter was defined as the spectrum of the total 
magnetic fluctuation variance {SB\). V{oj), however, is de- 
fined as the spectrum of the total fluctuation energy density 
(see eqs. IfTSl and Il22ll below). 

Cranmer et al. (1999a) assumed that the power spectrum 
scaled with the overall radial dependence of {SB\), as com- 
puted via MHD wave action conservation. However, the rel- 
atively small amount of energy in the high-frequency tail of 
the spectrum may not necessarily scale with the total inte- 
grated power. Thus, here we modify this radial dependence 
using the MHD turbulence theory of Cranmer & van Balle- 
gooijen (2003) to take account of an assumed anisotropic cas- 
cade that maintains the high-frequency tail. We assume that 
7^(17,, r) is a separable function of radius and frequency (for 
which uj ri,). The radial dependence is given by 



P(a.,r) cx 



3/2 



P^/Hup + Va)' 



(10) 



where p is the mass density. This scaling relation is derived 
in Section |7]below, and the normalization constant to be used 
is the main adjustable parameter that we modify to produce 
agreement with the measured values of T±i at 1.069 Rq. 

The background plasma properties are given by the empir- 
ical model of Cranmer & van Ballegooijen (2005). These 
properties include the polar flux-tube area A(r), the outflow 
speeds (for which we assume that the ions are flowing at the 
same speed as the protons; m, = Up), and the number densities. 
The primary plasma number density (rip w rtg) is assumed to 
have the same radial dependence as the ion number density 
«,. Note that the absolute value of the ion number density is 
never needed, only its radial dependence. 

The Cranmer & van Ballegooijen (2005) model was essen- 
tially a "cold plasma" treatment, so our description of the 
background proton temperature Tp is given by scaling the self- 
consistent polar coronal hole temperature of Cranmer et al. 
(2007) up and down by arbitrary constant factors. The loca- 
tion of the transition region was also scaled down to a stan- 
dard height of 0.003 Rq above the photosphere. We assume 
isotropic protons with Ty = T±p = Tp. The actual proton tem- 
perature in the corona is not well constrained by observations, 
since line ratio techniques are sensitive to the electron temper- 
ature, and the widths of H" lines are more strongly "contami- 
nated" by non-thermal motions than the lines of heavy ions. 

In order to determine what to use for the proton tempera- 
ture in the models. Figure |6] shows a range of measurements 
and model predictions for Tp, as well as the electron tem- 
perature Te (which is expected to be collisionally coupled to 
the protons at high enough densities) and the mean one-fluid 
temperature, which should be approximately proportional to 
(Tp + Te)/2. Off-limb measurements of T^ in polar coronal 
holes are shown from SUMER data described by Wilhelm 
(2006) and Landi (2008). Proton temperatures at large heights 
have been measured by UVCS/SOHO, and the plotted points 
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Fig. 6. — Assembled temperatures in polar coronal holes: Mean plasma 
temperatures from a semi-empirical VAL/FAL-type model {dashed line; 
Avrett & Loeser 2008) and from a turbulence-driven coronal heating model 
{solid line; Cranmer et al. 2007). Te from off-limb SUMER measurements 
made by Wilhelm (2006) {dark blue bars) and Landi (2008) {light blue 
bars). Tp from off-Umb SUMER measurements by Marsch et al. (2000) 
{circles), and UVCS measurements assembled by Cranmer (2004), with 
data from Cranmer et al. (1999b) {squares), Esser et al. (1999) {diamonds), 
ZangrilH et al. (1999) {asterisks), and Antonucci et al. (2000) {triangles). 

are those assembled by Cranmer (2004), with a model pre- 
diction for the non-thermal wave amplitudes subtracted from 
the empirical kinetic temperatures. The transition-region-like 
values of Tp at low heights were measured by Marsch et al. 
(2000); these are quite low compared to most time-steady 
coronal models that place the sharp transition region at a lower 
height around w 0.003 Rq. These observations may have 
been influenced by the prevalence of narrow and cool spicules 
above the limb and may not be representative of the coronal 
footpoints of solar wind streams. 

Before proceeding to discuss the solution of equation (|5]l, 
it is important to make clear the relative strengths of the in- 
dividual heating and cooling terms in this conservation equa- 
tion. As an illustration, Figure|2]shows the radial dependence 
of timescales computed from each term in equation (|5]l for the 
standard choice of parameters (and the O^^ ion) discussed in 
§ 8. Given that the terms in this equation have the units of 
a spatial temperature gradient, a timescale can be constructed 
by multiplying by the quantity Ui/T±i and taking the inverse. 
For example, the timescale corresponding to ion cyclotron 
heating is given by 



Theat 



nikBT_ 



Q±i 



(11) 



The timescale for adiabatic cooling Tadia is given by using 
the second term on the left-hand side of equation (jSj. Two 
timescales for Coulomb collisions can be constructed by us- 
ing either the ion temperature (i.e., the final term in eq. ||5|) or 
the proton temperature (the pentultimate term). These gener- 
ate the collisional timescales Tcoh.; and Tcou.p, respectively. As 
can be seen from Figure |7] at low heights the most rapid pro- 
cesses are the Coulomb collisions. Above a height of approx- 
imately 0.1 Rq, the ion cyclotron heating timescale begins to 
be as rapid as collisions, and the ion temperature T±i begins to 




1.00 



Fig. 7. — Radial dependence of timescales coiresponding to individual 
terms in the energy conservation equation (3); see text for definitions. 

increase above the proton temperature. At still larger heights 
in the extended corona (r w 2Rq) adiabatic cooling begins to 
be as important as ion cyclotron heating, and Coulomb colli- 
sions are unimportant. 

The numerical integration technique used to compute T±i{r) 
from equation (|5]l was slightly different from that used by 
Cranmer et al. (1999a). Because of the strong collisions near 
the coronal base, the earlier use of simple Euler steps became 
highly unstable for realistic grid spacings. Instead, we used a 
piecewise analytic solution of the energy equation that takes 
account of the collisional coupling explicitly. Equation (|5]l 
can be expressed schematically as 



dT±i 
dr 



+ pT±i = q 



(12) 



where p and q are considered to be slowly varying functions 
of radius. The Coulomb collision term has been split up be- 
tween p and q. The analytic solution of this equation, which 
we use to step the solution for T±i up from height rj to height 
rj+i, is given by the integrating factor method as 



i±i,j+i - 



= — + e 



-p{rj+i-rj) 



(13) 



In the limit of weak collisions (/? <C 1), this reduces to the 
standard Euler step. 



T±i,j + q{rj+i-rj) . 



(14) 



In the limit of strong collisions (i.e., ignoring the adiabatic ex- 
pansion and cyclotron heating terms), the ratio q/p reduces to 
the proton temperature Tp, and equation ( fT3] ) is seen to stably 
drive the solution towards thermal equilibrium. 

7. EXPECTATIONS FROM ANISOTROPIC TURBULENCE THEORY 

As an independent prediction of the resonant wave power, 
we aim to describe the anisotropic evolution of wave energy 
in 2D wavenumber space (k\\,k±) by means of the advection- 
diffusion theory described by Cranmer & van Ballegooijen 
(2003) and in Appendix C of Cranmer & van Ballegooijen 
(2005). The theoretical development given in this section is 
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somewhat independent of the "semi-empirical" derivation of 
wave power that comes from solving the equations of the pre- 
vious section. 

The goal is to develop an analytic solution for the turbulent 
power spectrum (for frequencies near the ion cyclotron res- 
onances) in the limit that the cascade is allowed to proceed 
to its final "driven" steady state in a small homogeneous vol- 
ume of plasma. The three-dimensional total power spectrum 
is defined as 



Ml 

47r 



(15) 



where f/tot is the total energy density of the Alfven waves, and 
we write the volume element d^k in cylindrical coordinates 
as 2TTk±dk±dk^^ to assume symmetry in the two directions 
perpendicular to the background field. 

For fully developed anisotropic MHD turbulence, we as- 
sume the power spectrum to be a separable function of two 
variables: k± and a nonlinearity parameter y defined as the 
ratio of the local wind-frame frequency VaA;|| to an assumed 
nonlinear eddy turnover rate {SV)k± (see, e.g., Goldreich & 
Sridhar 1995). We use the notation from Section 2.3 of Cran- 
mer & van Ballegooijen (2003) and define 



P3D(k\\,k±) 



g(y) 



(16) 



where W±(k±) is a reduced power spectrum that describes the 
dominant perpendicular cascade. For the MHD inertial range, 
the reduced spectrum is given as 

W (h. \ - j Utot(kj_/koM)~^^^(3Trp)-^ , kout <k±< kn 
yyjLKi^n <^Q^ otherwise 

(17) 

The above definition assumes the existence of a finite outer- 
scale perpendicular wavenumber ^out, which we assumed to 
be inversely proportional to the correlation length of the tur- 
bulence, and an inner-scale wavenumber kin ^ ^out, which we 
assume to be equivalent to the inverse proton gyrofrequency. 
The factor of 3ttp above is needed to normalize the full power 
spectrum as defined in eq. (flST l. 

The A: 1 1 dependence of the power spectrum is contained in 
the dimensionless g{y) function given in eq. ( fTSI l. We define 



y = 



k\\VA 



k±W^ 



1/2 



(18) 



The condition y = 1 is defined as "critical balance" by Goldre- 
ich & Sridhar (1995), and their analysis only constrains the 
general shape of g(y), not its exact value. The critical balance 
condition captures the highly nonlinear state of turbulence, for 
which a coherent wave survives for no more than about one or 
two periods before nonlinear processes transfer its energy to 
the smaller scales. Cranmer & van Ballegooijen (2003) solved 
a simple wavenumber diffusion equation to obtain an analytic 
relation for g(y). For the MHD inertial range, this relation can 
be shown to be equivalent to 



g(y) = 



2r(«) 



3r(n-0.5)V^ 



1 + 



4/ 



(19) 



which is normaUzed to unity when integrated over all y, and 
with 

33 

«=! + /. (20) 
47 



The dimensionless constants (3 and 7 describe the relative 
strengths of advection and diffusion, respectively, in the k± 
direction. They only occur as the ratio /3/7, which we take 
to be a free parameter. Cranmer & van Ballegooijen (2003) 
discussed the most realistic values for this ratio, which is rel- 
atively unconstrained by existing turbulence simulations. The 
"random walk" turbulence model of van Ballegooijen (1986), 
though, suggested that /3/7 ~ 1 . However, Cranmer & van 
Ballegooijen (2003) found that one would need this ratio to 
be smaller than about 0.25 in order to produce enough paral- 
lel cascade in the corona to heat protons and heavy ions via 
cyclotron resonance. 

With the above definitions it becomes possible to integrate 
over k± and obtain a reduced ID power spectrum 



Pid(A||) = 2tt I dk^k^Piu{kn,k±) 



(21) 



which is related to the frequency spectrum Viui) via 

PiD(k\\)dki\ = V(uj)dLu . (22) 

If we assume, for simplicity, that the ion cyclotron waves of 
interest are far enough from the proton resonance for the stan- 
dard MHD dispersion relation (uj = A:||Va) to hold, the above 
parameterizations for ^30(^11 ,A:^) can be applied to show that 

giy) 



47rpf/tot, 1/3 



dk I 



(23) 



This equation, combined with the cyclotron resonant fre- 
quency condition 

u}-Uik\\-Vli = Q , (24) 

is used to produce the "theoretical" power spectra (as a func- 
tion of the /3/7 parameter) that are compared with the empir- 
ically derived power spectra at 1.069 /?©. For ion cyclotron 
waves at any given height in the low corona (where m, <C Va), 
the wave frequency uj is directly proportional to the ion charge 
to mass ratio Zi/Aj. The main plasma quantities, including the 
Alfven wave scale lengths fcout and k\a, the total energy den- 
sity f/tot, and the ion Larmor frequencies il,, are taken from 
the radially dependent coronal hole model of Cranmer & van 
Ballegooijen (2005). 

The above equations specify how the resonant frequency 
spectrum V{uf) can be computed, but it is also helpful to de- 
rive a concise relation for the large-scale radial dependence of 
the wave power in the ion cyclotron frequency range. The 
gyroresonant wave power exists in the limit of large val- 
ues of the dimensionless critical-balance parameter y. For 
y ^ 1, the function g(y) cx y~^", and this implies that the high- 
frequency tail of the spectrum has a similar power-law depen- 
dence V oc bj~^". Cranmer et al. (1999) treated this exponent 
as a free parameter 77. Utilizing the definitions given above, 
this proportionality can be expressed in terms of more basic 
plasma parameters as 



giy) oc li 



4«/3 



( ^out^tot A 
3^p02 j 



(25) 



Note that the actual frequency dependence is "hidden" in in 
the ion gyrofrequency fi,, which is simply Z,/A, times the 
proton gyrofrequency 51^. Performing the integration in equa- 
tion ( I23I 1 and removing all constants thus provides the radial 
dependence of the power spectrum: 

„0.5-«jy-«+0.5;,(2n+l)/3 



V cx 



2(n+2)/3 4(«-l)/3 



(26) 
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Fig. 



Radial dependence of modeled proton temperature (thin solid 
line) and perpendicular O^^ temperatures, assuming a range of resonant wave 
power normalization values given relative to Pq, the iterated value required to 
heat the ions to the observed level at 1.069 Rq (thick solid line). Resulting 
values for T±j(r) for power levels of IOPq (thin dotted line), 2Pq (dot-dashed 
line), 0.5Po (dashed line), and O.lPo (thick dotted line) are also shown. 

where the proton thermal speed Wp is part of the definition 
of the mean proton gyroradius used in ki-a- Also, the constant 
ratio of Zj/A,- has been removed in order to isolate the radial 
dependence for any given ion resonance. 

We will see below that the limiting case of /3/7 <C 1 (i.e., 
n sa 1, or ry sa 2 in the notation of Cranmer et al. 1999) may 
be appropriate for comparison with the observations, so let us 
apply this limit to obtain 



(27) 



The outer scale wavenumber ^out is usually assumed to be the 
reciprocal of the perpendicular correlation length, and the lat- 
ter is often taken to be proportional to the flux-tube cross sec- 
tion, so we assume A:out oc A"'/^ cx b\I^ (Hollweg 1986). We 
also know that Hp oc Bq. Finally, the total Alfven wave energy 
density can be assumed to scale with the WKB conservation 
of wave action, which is given by 



f/tc 



Va 



A{up + VaY 



(28) 



Applying these scaUngs then leads to the radial dependence 
for the resonant power given in equation ( fTOl i. 

8. MODEL RESULTS 

The semi-empirical equation of ion energy conservation 
(eq. 15] ) was integrated numerically from the coronal base 
(1.003 Rq) to the reference radius of 1.069 Rq for each of 
25 ions having measured SUMER line widths at this height. 
These ions have Z,/Ai values ranging between 0.161 (Fe^""") 
and 0.370 (Mg^"^), and ion temperatures between 1.69 MK 
(Ne*""^) and 7.22 MK (Fe"+). Also, several grids of models 
were created for a range of proton temperatures, all with the 
same relative radial dependence but with values at 1 .069 Rq 
ranging between 0.5 and 3 MK. 
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• 0.75 MK 

• 1.00 MK 

• 1.25 MK 

• 1.50 MK 

• 1.75 MK 
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Fig. 9. — Ion cyclotron wave power at 1.069 Rq plotted versus ion charge- 
to-mass ratio Z,/A, (the latter a proxy for frequency). Empirically determined 
wave power levels for individual ions (filled circles; see color key for Tp 
normalization) are compared with theoretical predictions for the wave power, 
based on adopted radially dependent models of the total wave power, Alfven 
speed, and other plasma properties (black solid cun'es). 

Figure [8] illustrates the result of this numerical integration 
for a representative ion (O^"*") and a specific choice for the 
proton temperature. The optimal (iterated) wave power nec- 
essary to heat these ions to the observed temperature at the 
reference height is denoted Pq. Figure[8]also shows the result- 
ing T±i{r) curves that are obtained when Pq is varied up and 
down by factors of 2 to 10. The ion temperature can become 
several orders of magnitude larger than the proton tempera- 
ture at larger heights. This has been confirmed by the UVCS 
O VI measurements between 1.5 and 3.5 Rq (e.g.. Kohl et al. 
2006; Cranmer et al. 2008). Note, though, that at such lai-ge 
heights, where Coulomb collisions become very weak, a more 
proper treatment of the resonant wave-particle interactions 
is warranted — either via the coupled T±i and 7]|, equations 
(Cranmer et al. 1999a) or even a fully kinetic treatment of the 
ion velocity distributions (e.g., Isenberg & Vasquez 2007). 

Figure |9] shows the result of repeating the numerical inte- 
gration and wave-power iteration for all 25 ions at 1.069 Rq. 
Recall that no assumptions were made about the frequency 
dependence of the resonant wave power spectrum for these 
models. The power for each ion was computed independently 
of the others. The resonant wave power is shown as a function 
of Zj/A, for each of the 16 choices for the Tp normalization. 
Straight-line fits are shown for two distinct regions of Z,/A, 
values, in order to guide the eye and make rough comparisons 
to the theoretical expectations for Viuj). The black curves 
in Figure |9] are shown to compare the semi-empirical results 
with the expectations of turbulence theory. In other words, 
these curves come from evaluating equation ( |23] ) numerically 
for V{ijj), assuming a range of j3/^ values between 0.01 and 
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1. 

For large enough choices for Tp, there occur several ions 
that have smaller perpendicular temperatures than the adopted 
proton temperature at 1 .069 Rq . These are shown with dashed 
lines for the straight-line fits, and thus should not be consid- 
ered as self-consistent models for the entire ion data set. In 
Figure |9] the "green" curves indicate a collection of models 
all with Tp^ 2.1 MK, because we wanted to explore specif- 
ically the behavior of the O^"^ heating in the limit of strong 
collisionality (i.e., because « 2.15 MK for this ion); see 
below. 

For the lowest frequencies (i.e., Z,/A, < 0.27), both the 
spectral slope and the absolute wave power at 1.069 Rq seem 
to agree well with the turbulent cascade predictions for val- 
ues of P/'j ^ 0.3. This agreement is actually quite notable, 
since the theoretical predictions (black curves) were normal- 
ized by pre-existing assumptions about, e.g., the total wave 
power f/tot, and the empirical sets of points were iterated com- 
pletely independently from any considerations of MHD turbu- 
lence. If the proton temperature in the low corona remains less 
than about 2 MK (see, e.g., Figure|6]l, then this result provides 
some additional observational justification for a value of f3/j 
of order 0.1 to 0.3 (see Cranmer & van Ballegooijen 2003). 

For the highest frequencies {Zi/Ai > 0.27), it seems clear 
that some kind of "upswing" in the resonant wave power is 
indicated by the empirical points in Figure |9] If this is indeed 
a real effect, there are several possible explanations in terms 
of MHD and kinetic fluctuations: 

1. The additional power at large gyrofrequencies may 
arise because of plasma instabilities centered around 
the cyclotron resonances of alpha particles {Zi/Ai = 0.5) 
or protons (Z,/A,- = 1); see, e.g., Markovskii (2001), 
Isenberg (2001), Zhang (2003), Laming (2004), and 
Markovskii et al. (2006). Thus, the effects seen at 
charge-to-mass ratios around 0.35 may just be the tail 
of an additional population of waves that peak at larger 
frequencies. 

2. This may be a similar effect as the spectral flattening 
(i.e., enhancement above a power law frequency de- 
pendence) evident in radio scintillation measurements 
of density fluctuations at larger heights (e.g.. Coles & 
Harmon 1989). This flattening was modeled success- 
fully by Harmon & Coles (2005) by taking into account 
the enhanced compressibility of obliquely propagating 
Alfven waves once they reach ion cyclotron frequencies. 
It may be possible for compressibility effects to either 
alter some aspects of the turbulent cascade or change 
the functional dependence of the heating rate on the 
available wave power (i.e., eq. lID) in order to produce 
this effect. 

3. A somewhat more speculative idea is that the particu- 
lar type of MHD turbulence acting in the low corona 
may undergo a kind of bottleneck effect where the wave 
power piles up near the dissipation range. This effect 
appears in many turbulence simulations, and it has of- 
ten been suspected of being numerical in origin (e.g., 
Verma & Donzis 2007). However, there have been in- 
dependent theoretical predictions of similar kinds of 
wave-power pileup arising from nonlocal interactions 
between disparate scales in ^-space (e.g., Falkovich 
1994; Biskamp et al. 1998). 
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Fig. 10. — Radial dependence of modeled proton temperatures {dotted 
lines) and pemendicular O^^ temperatures (solid lines), with the same color 
key as Figure|9] Also shown are observational lower and upper limits on T^i 
from SUMER and UVCS (black solid lines). 

Figure [To] illustrates the radial dependence of the ion heat- 
ing for O^"^ ions, and shows how 7j^, varies with the different 
choices for Tp. The color scheme for proton temperature val- 
ues is the same as in Figure|9] The black curves at the heights 
corresponding to the SUMER measurements show the lower 
and upper limits defined by equation The black curves 
at larger heights give the limits on the ion temperature from 
recent UVCS O VI empirical models (Cranmer et al. 2008). 

The models with relatively low proton temperatures {Tp<2 
MK) require a large amount of resonant wave power to "com- 
bat" Coulomb collisions and raise the ion temperature to the 
observed SUMER value at 1.069 Rq. Unfortunately, this 
large resonant wave power then causes the ion temperature 
to "explode" to unrealistically huge values at larger heights 
where the plasma becomes collisionless. On the other hand, 
the green curves in Figure [TO] agree reasonably well with both 
the SUMER and UVCS data and correspond to larger as- 
sumed values for Tp. Here, the reason for a fine grid of proton 
temperature curves having 

0<^-l«l (29) 

becomes evident. These choices for Tp do not need as much 
wave power to reach the observed SUMER ion temperature 
at 1.069 Rq, and thus they are less likely to result in unre- 
alistically large temperatures at larger heights. However, it 
is important to note that these best-fitting solutions seem to 
have quite strong Coulomb collisions even at the lowest of 
the SUMER heights. This would imply that at, say, 1.03- 
1.05 Rq, all of the ion temperatures would be expected to be 
tightly coupled to both Tp and to one another This, as seen in 
Figure |3] above, is not the case. Clearly the models presented 

^ It is possible, of course, that Coulomb collisions are not the only means 
of thermalizing and isotropizing the ions. The existence of various kinds of 
turbulent fluctuations has been suggested to be able to can couple particles 
and possibly provide extra "quasi-collisions" (e.g., Perkins 1973; Dum 1983; 
Kellogg 2000). This may need to be investigated in the context of the low 
corona. 
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here are only a first step toward a self-consistent description 
of the preferential ion heating that applies for the full range of 
heights observed by SUMER and UVCS. 

9. CONCLUSIONS 

In the present work we have used SUMER observations of 
a polar coronal hole to measure the ion temperatures 7] for a 
large number of ions, in order to determine their dependence 
on the charge-over-mass ratio Zj/A,. We repeated our mea- 
surements for several heights between 1.03 Rq and 1.17 Rq 
in order to investigate the height dependence of our results. 
We used the method devised by Tu et al. (1998) to determine 
Ti; this method only assumed that non- thermal velocities v„,h 
are the same for all ions. 

We found that ion temperatures are larger than the electron 
temperatures in nearly all cases and at all heights, and that 
results are qualitatively the same in the 1.03-1.17 Rq range 
of distances from the hmb, although the number of available 
ions decreases as distance from the limb increases. Our most 
notable result is an U-shaped dependence of 7] from Z,/A„ 
where ions with low Z,/A, (less than about 0.23) have very 
large T; values, while those with Z,/A,- > 0.23 have roughly 
constant 7). 

We used our measured Ti values to constrain an exploratory 
model of solar wind heating and acceleration based on ion- 



cyclotron waves, under the assumption that such waves are 
gradually replenished as they are dissipated. The measured 
Ti values help us constrain the proton temperatures as well as 
the ratio between advection and diffusion (in a likely scenario 
of anisotropic MHD turbulence). We find that the observa- 
tions of ions having Z, /A, values smaller than about 0.25 are 
consistent with a turbulence model very similar to that shown 
by Cranmer & van Ballegooijen (2003) to be able to energize 
protons in the extended corona. 

However, 7; measurements obtained for ions with Z,/A,- 
larger than 0.25 show an upswing in wave power that is 
difficult to reconcile with traditional views of turbulent cas- 
cade. We discuss the implications of this inferred increase in 
wave power, suggesting several different possibilities such as 
plasma instabilities, obliquely propagating Alfven-wave com- 
pressibility effects, or a turbulent bottleneck effect. 
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like to thank Adriaan van Ballegooijen and Mari Paz Miralles 
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